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Abstract 



Annihilations of weakly interacting dark matter particles provide an 
important signature for the possibility of indirect detection of dark mat- 
^ ' ter in galaxy halos. These self-annihilations can be greatly enhanced in 

the vicinity of a massive black hole. We show that the massive black hole 
present at the centre of our galaxy accretes dark matter particles, creating 
a region of very high particle density. Consequently the annihilation rate 
is considerably increased, with a large number of e + e~ pairs being pro- 
duced either directly or by successive decays of mesons. We evaluate the 
, synchrotron emission (and self-absorption) associated with the propaga- 

— ■ tion of these particles through the galactic magnetic field, and are able to 

constrain the allowed values of masses and cross sections of dark matter 
particles. 
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to 1 Introduction 

There is convincing evidence for the existence of an unseen non-baryonic com- 
ponent in the energy-density of the universe. The most promising dark matter 
candidates appear to be weakly interacting massive particles (WIMPs) and in 
particular the so-called neutralinos, arising in supersymmetric scenarios as well 
as much heavier particles (WIMPzillas) which could have been produced non- 
thermally in the early universe (for a review of particle candidates for dark 
matter see e.g. Ellis (1998)). 

The annihilations of these X-particles would generate quarks, leptons, gauge 
and Higgs bosons and gluons. Consequently e + e~ pairs would be produced ei- 
ther directly or by successive decays of mesons, and they are expected to lose 
their energy through synchrotron radiation as they propagate in the galactic 
magnetic field. This radiation is expected to be greatly enhanced in the prox- 
imity of the galactic centre, where the existence of a massive black hole creates 
a region of very high dark matter particle density and consequently a great 
increase in the annihilation rate and synchrotron radiation. 

We discuss in section 2 the distribution of dark matter particles in our galaxy 
and in particular around the central black hole, following Gondolo & Silk (1999) 
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(from now on, Paper I). In section 3 we evaluate the annihilation e + e _ spectrum 
and their synchrotron emission. Finally in section 4 we discuss our results and 
determine the regions of the mass-annihilation cross section plane which give 
predictions compatible with experimental data. 



2 Dark matter distribution around the galactic 
centre 

There is strong evidence for the existence of a massive compact object lying 
within the inner 0.015pc of the galactic centre (see Yusef-Zadeh, Melia & Wardle 
(2000) and references therein). This object is a compelling candidate for a 
massive black hole, with mass M = (2.6 ± 0.2) x 10 6 M Q . The galactic halo 
density profile is modified in the neighborood of the galactic centre from the 
adiabatic process of accretion towards the central black hole. If we consider an 
initially power-law type profile of index 7, similar to those predicted by high 
resolution N-body simulations (Navarro, Frenk & White 1997; Ghigna et al. 
2000), the corresponding dark matter profile after accretion is, from paper I 



a. 
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where 7s P = (9— 2j)/(4— 7), D is the solar distance from the galactic centre and 
Pd = 0.2AGeV/ c 2 /cm 3 is the density in the solar neighbourhood. The factors a 7 
and g~f(r) cannot be determined analytically (for approximate expressions and 
numerical values see paper I). The expression (|lj) is only valid in a central region 
of size R sp = a 1 D{M / prjD 3 ) 1 ' ( 3 ^ where the central black hole dominates the 
gravitational potential. 

If we take into account the annihilation of dark matter particles, the density 
cannot grow to arbitrarily high values, the maximal density being fixed by the 
value 

Pcore = m/avtBH (2) 

where tBH ~ 1Q W years is the age of the central black hole. The final profile, 
resulting from the adiabatic accretion of annihilating dark matter on a massive 
black hole is 

I \ P'(r)pcore 

Pdm{r) = , (3) 

P(r) + Pcore 

following a power law for large values of r, and with a flat core of density p core 
and dimension 

Rm = R J^\ l ' M (4 ) 
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Figure 1: Values of Y e as a function of particle mass. 

3 Constraints from Synchrotron Emission 

Among the products of annihilation of dark matter particles, there will be elec- 
trons and positrons, which are expected to produce synchrotron radiation in the 
magnetic field around the galactic centre. The e + e~ component produced by 
the hadronic jets has been computed using the MLLA approximation^]. 

The galactic magnetic field can be described using the 'equipartition as- 
sumption', where the magnetic, kinetic and gravitational energy of the matter 
accreting on the central black hole are in approximate equipartition (see Melia 
(1992)). In this case the magnetic field can be expressed as 

B(r) = 1/iG -j (5) 



Energy-loss length scale. As we shall see, most of the annihilations happen 
at very small distances from the centre, typically at w min(R core ,lORs), i.e. 
in a region with magnetic fields of the order of > 1G. Under these conditions, 
comparable to the size of the region where most of the annihilations occur, the 
electrons lose their energy almost in place. To see this, consider the critical 
synchrotron frequency v c (E~) i.e. the frequency around which the synchrotron 

iFor details see Dokshitzer et al. (1991), Ellis, Stirling, & Webber (1996), Khoze k, Ochs 
(1997). For applications to ultra-high energy cosmic rays see Bhattacharjee & Sigl (2000). 
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Figure 2: Left panel: A v as a function of frequency for mx — ITeV. The two 
upper curves correspond to the cross section av w 10~ 28 /m^-(GeV) cm 3 s _1 , 
close to the unitarity limit; the two lower curves correspond to av m 
10~ 38 /m^(GeV) cm 3 s _1 , a cross section more typical for wimps. Results for 
two values of the density profile are shown in each case: 7=1 (solid curves) 
and 7 = 1.5 (dashed curves). Right panel: A v as a function of the particle mass 
for z/=408MHz, av = 10 -10 /m x (in physical units) and two values of 7. 



emission of an electron of energy E, in a magnetic field of strength B, peaks, 
namely 

3 eB ( E \ 2 

v c {E) = - 2 . (6) 

47T m e c \m e c J 

Inverting this relation, we find the energy of the electrons which give the maxi- 
mum contribution at that frequency, 

The typical synchrotron loss length for an energy corresponding to E m {v) can 
be expressed as a function of the distance r from the central black hole for the 
magnetic field profile in Eq. (|5|): 



V2 /„5J1\ 1/2 

l e {E m { V )) 



( 27 y /m a e c L1 Y'" 1 1 
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For a frequency of 408 MHz, which produces, as we shall see, the most stringent 
upper bound on the Galactic Centre emission, we get 

/ \ is/8 

le{E m {y)) = 5.074 • 10 8 ( — J pc (9) 

The diffusion length D = D(B(r), E), which can be approximated as a third of 
the radius of gyration of the electron, is given by 



E fi / v r \ i/ 2 / r 



15/8 

^M) = ^«8.991.10-^j- {-) pc (10) 

Since l e (E m (v)) » D{E m {i>)) for all practically relevant parameters, an elec- 
tron will diffuse a distance \Jl e D before it loses most of its energy. Numerically, 



w 3.035- 10 1 — 
r \pcj 



7/8 

(11) 



or, expressing the distance as a function of the Schwarzschild radius 



CD 



r 



7/8 



7.059 • 10~ 5 ( — ) (12) 



We can thus assume that the electrons lose their energy practically in place. 

Electron Production Spectrum. To compute the synchrotron luminosity 
produced by the propagation of e^- in the galactic magnetic field, we need to 
evaluate their energy distribution in the magnetic field which is given by (see 
Gondolo (2000)) 

dn _ TY e (> E) 



dE P(E) 
where Y is the annihilation rate 



-fe(r) (13) 



/•OG 

r = 5/ dr. (14) 

m x JO 

The function f e (r) is given by 

= f°° f? 2 , (15) 

Jo Ps P ^ r dr 

and 
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Figure 3: Comparison of Sgr A* (see Narayan et al. (1998)) observed spectrum 
with expected fluxes. The values of particle mass and cross section were choosen 
to fit the experimental data normalisation. 



is the total synchrotron power spectrum. Note that the general expression for 
f e {r) would have to take into account spatial redistribution by diffusion (see 
e.g. Gondolo (2000)), which we demonstrated to be negligible in our model. 

The quantity Y e (> E) is the number of e + e~ with energy above E produced 
per annihilation, which depends on the annihilation modes. Eq. (]?]) shows that 
for the frequencies we are interested in, E m (y) « mj, and thus the energy 
dependence of Y e (> E) can be neglected. We will estimate Y e (> E) by the 
number of charged particles produced in quark fragmentation (see footnote) ; in 
figure [I] we show the values of Y as a function of the particle mass. 

Synchrotron Luminosity. For each electron the total power radiated in 
the frequency interval between v and v + dv is given by 

./q p 3 /•oo 



V3p 3 



where we introduced 

F 



v c {E) 



G 



Integrating this formula we obtain the total synchrotron luminosity 

dr Avr 2 / dE ^P{v,E) (19) 
Jm e d E 

which by substitution becomes 



It is possible to simplify this formula by introducing the following approxi- 
mation for the function F (see Rybicki & Lightman (1979)) 



F [ ^ ) « 5{u/u c {E) - 0.29) . (21) 



The evaluation of the integral then gives 



8 V0.297T e 

where 



9 / 1 m 3 p c 5 \ 1/2 TYJ> E) , N 

L - ~ o 7^ > ^ / (22) 



ATir 2 f e {r)B- 1/2 {r) (23) 



Synchrotron Self- Absorption. The synchrotron self-absorption coefficent 
is by definition (see Rybicki & Lightman (1979)) 

1 f°° 

A v = — (1 - e _T W)27r& db (24) 
where r(6) is the optical depth as a function of the cylindrical coordinate b 

/>oo 

7-(&) =a v / / e (M) (25) 

Jd(b) 



and the coefficent o„ is given by 



e 3 TB(r) f m E 2± \Ye(> E) 



9m e is 2 /_ <iE 



E 2 P{E) 



F I — ) rf£ (26) 



The final luminosity is obtained by multiplying Eq. ( |19[ ) with A„ given by 
Eq. (24). It is evident that in the limit of small optical depths the coefficent 
A v — > 1, as can be seen by expanding the exponential. 

Furthermore the lower limit of integration of expression |25| is 

d(b) = for b 2 +z 2 > (4Rs) 2 (27) 



d(b) = v (ARs) 2 — 6 2 elsewhere 



7 




10° 10 2 10 4 10 6 10' 

Mx (GeV) 



Figure 4: Exclusion plot based on the comparison between predicted flux and 
radio observations of the galactic centre. The 3 solid curves indicate, for different 
values of the density profile power law index, the lower edge of the excluded 
regions. The dashed line shows, for comparison, the unitarity bound, av ~ 
\/m 2 x . The shaded region is the portion of the parameter space occupied by 
cosmologically interesting neutralinos (i.e. those leading to 0.025 < ilxh 2 < 1; 
see, e.g. Bergstrom, Ullio & Buckley (1997)). 



Using the approximation introduced in Eq. (|21j) we find for a v the following 
expression 



TY - 2 



, : - (28) 

which can in turn be used to evaluate r(b) in Eq. ( jga ) and A v in Eq. (pij). 

Note that the synchrotron loss time-scale is proportional to B~ 2 oc r^> 2 . We 
compare this with the annihilation time, which determines where the inner den- 
sity drops drastically: t ann oc r' 9_27 ^( 4_7 ^ +1 / 2 . In other words, the annihilation 
time goes to zero more rapidly than the synchrotron loss time as r goes to 0. 
Hence synchrotron losses are important throughout the annihilation region. 

Extension to the wimpzilla mass range introduces far more uncertain physics. 
The wimpzilla cross sections strongly depend on new physics beyond the elec- 
troweak scale. Similarly to the case of electroweak gauge bosons, one can expect 
a general scaling av ~ a 2 /m 2 x for the self-annihilation cross section where a < 1 
is some dimensionless gauge coupling. In four dimensional field theory, the con- 
servation of probability (unitarity) roughly corresponds to this scaling for a — * 1 
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(see, e.g., Weinberg (1995)). In order to retain generality, we will therefore con- 
sider Wimpzilla annihilation cross sections in the range av < l/m 2 x . We note, 
however, that beyond four dimensional field theory, such as in theories with 
extra dimensions and in string theory, larger cross sections may be possible. 

We first evaluate the self-absorption coefficent for selected values of the mass, 
as a function of frequency. The result, shown in figj^, is a coefficent which grows 
from very low values (showing that absorption is important at 408 MHz) and 
then reaches the value 1, around a frequency which is strongly dependent on the 
cross section and the mass mx , but not so much on the profile power law index 
7. On the left the above coefficent is evaluated for two different values of the 
cross section, the first one corresponding to the maximum possible value (the so 
called unitary bound, see Griest & Kamionkowski (1990)) and the second one 
corresponding to typical cross sections one can find in supersymmetric scenarios 
(see, e.g., Bergstrom, Ullio & Buckley (1997)). 

The right part of figj| shows the self-absorption coefficent at the fixed fre- 
quency of 408MHz as a function of the neutralino mass. The behaviour shown 
is qualitatively the same for any value of the cross-section and for different 7. 
In figure ^ we compare the predicted spectrum with the observations; we choose 
a set of parameters mx, 7 and av in order to reproduce the observed normal- 
isation. It is remarkable that in that way one can reproduce the observe d 
spectrum over a significant range of frequencies. The set of dark matter param- 
eters for which the fluxes predicted in our model are consistent with observation 
is shown in the exclusion plot of fig. |^. The boundaries of the excluded range 
represent the parameter values for which the observed flux is explained by the 
dark matter scenario discussed here. 



4 Conclusions and perspectives 

We have shown that present data on the emission from Sgr A* are compatible 
with quite a wide set of dark matter parameters. The evaluation of synchrotron 
self-absorption has enabled us to reach an alternative conclusion from an ear- 
lier claim of incompatibility of cuspy halos with the existence of annihilating 
wimp dark matter. Our final result is that the experimental data on Sgr A* 
spectrum at radio wavelengths could be explained by synchrotron emission of 
electrons produced in the annihilation of relatively massive dark matter parti- 
cles, extending from TeV masses to mx > 10 8 GeV. The former is relevant to 
recent studies of coannihilation (Boehm, Djouadi & Drees (2000)), that suggest 
that WIMPs with Qxh 2 = 0.2 can extend up to several TeV, and the latter 
is relevant for particles (wimpzillas) that are produced non-thermally in the 
primordial universe. 

We have found that with the current data situation, the synchrotron emis- 
sion tends to give somewhat sharper constraints on masses and cross sections 
than the gamma-ray fluxes (cf. Baltz et al. (2000)). This situation could be 
reversed by more sensitive gamma-ray observations anticipated from upcom- 
ing experiments. However the synchrotron predictions are uncertain because of 
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our relative ignorance about the magnetic field strength near the central black 
hole, and gamma ray fluxes are subject to similarly uncertain amounts of self- 
absorption. The ANTARES neutrino experiment will eventually be able to set 
a relatively model-independent limit on the annihilation flux from the Galactic 
Centre. 
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